Studying the anisotropy of the gravitational wave stochastic background with LISA 
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A plethora of gravitational wave stochastic backgrounds populate the sensitivity window of the Laser 
Interferometer Space Antenna. We show that LISA can detect the anisotropy of the background 
corresponding to the multipole moments of order I — 2 and 4. The signal-to-noise ratio generated by 
galactic white dwarf binary systems could be as high as ~ 60 for 3 yrs of integration, and LISA could 
provide valuable information on the spatial distribution of a variety of galactic sources. We also 
show that the cross-correlation of the data sets from two interferometers could marginally lead to 
meaningful upper-limits on the degree of isotropy of the primordial gravitational wave background. 
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Gravitational wave (GW) stochastic backgrounds gen- 
erated in the early-Universe and produced by large popu- 
lations of faint astrophysical sources are among the most 
interesting signals for GW laser interferometers. Their 
detection will provide a new arena for early-Universe cos- 
mology and compact-object astronomy. The Laser Inter- 
ferometer Space Antenna (LISA) represents the key 
instrument in the low-frequency band of the spectrum 
10^^ Hz — 10^^ Hz), complementary to Earth-based 
detectors, operating in the window 1 Hz — 10 kHz. 
LISA could prove to be the most suitable device to search 
for, and study stochastic signals: (i) Several cosmological 
models predict a primordial GW stochastic background 
with a flat energy density spectrum Qg^{f) ~ const 
(see 1^, and references therein, for a thorough review): 
indeed the spectral density of the signal Sh{f) generated 
at the detector output steeply increases at low frequen- 
cies as Sh{f) oc rigw(/)//'^; (ii) Stochastic backgrounds 
generated by the incoherent super-position of GWs from 
a variety of short-period binary systems (including close 
white dwarf binaries, W UMa binaries, and neutron star 
binaries) both galactic ||-||] and extra-galactic j^,^ are 
guaranteed signals for LISA, with several components 
that dominate the instrumental noise in the frequency 
range 10"'' Hz < / < 3 x lO^^ Hz. 

The key issue in searching for stochastic backgrounds, 
and gaining insights into their production mechanism and 
the underlying physics, is to identify unambiguously the 
GW signal in one single LISA data stream. In fact, a 
stochastic background is a random process which is in- 
trinsically indistinguishable from the detector noise. We 
have discussed this issue in Q, in the context of per- 
fectly isotropic signals, showing the advantage of cross- 
correlation experiments involving two data sets with un- 
correlated noise, and pointing out fundamental sensitiv- 
ity limits. More recently it has been shown that, by es- 
timating the LISA instrumental noise using the "sym- 
metrized Sagnac observable" [0 it is possible to de- 
tect an isotropic stochastic background with a sensitiv- 
ity close to the one that can be achieved through cross- 



correlation [ |ll| . However, several stochastic backgrounds 
are expected to be anisotropic. In particular, the back- 
ground generated by galactic sources carries a strong 
signature due both to the peripheral location of the so- 
lar system in the galaxy, and the spatial distribution of 
the sources. The primordial GW background, although 
expected to be intrinsically isotropic to a high degree, 
presents a clear dipole structure due to the proper mo- 
tion of the galaxy with respect to the cosmological rest 
frame. 

In this paper we show that the peculiar orbital con- 
figuration of LISA - the instrument baricenter follows 
an Heliocentric orbit with period T = 1 yr, and the de- 
tector plane, tilted by 60° with respect to the Ecliptic, 
counter-rotates with the same period T [|l| - can be ex- 
ploited to design a suitable data analysis scheme to detect 
and study anisotropic stochastic signals using the Michel- 
son configuration. This possibility was earlier suggested 
in p^ . Here we carry out a rigorous analysis from the 
point of view of data analysis, following and we ap- 
ply our results to the relevant astrophysical and cosmo- 
logical scenarios. We provide estimates for the optimal 
signal-to-noise ratio (SNR) as a function of the signal pa- 
rameters, and we show that, by detecting an anisotropic 
background, it is indeed possible to gain key insights on 
the spatial distribution of short-period binary systems. 
We also discuss the possibility of obtaining meaningful 
upper-limits on the degree of isotropy of the primordial 
GW background, and comment on the improvement of 
the sensitivity and information extraction that would 
come from the correlation of the data streams of two 
separated and suitably oriented instruments. 

We start by briefiy reviewing the fundamental data 
analysis concepts; we assume the reader familiar with 
the literature on the detection of stochastic backgrounds 
in GW observations |l3| - [l5[| . The GW background am- 
plitude at any time t and position x can be written in 
the transverse-traceless gauge as 
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^ab ('^' ^ ~ lj2,3) is the polarization tensor of the 
wave, hA(fi, f) are the two independent polarization 
amplitudes {A = +, x), and fl is the unit vector 
on the two-sphere along the wave's propagation di- 
rection. Here we consider stochastic signals that are 
Gaussian, unpolarized but not isotropic: {hA{f,^l)) = 

and {hA{f,n)h\,if',n')) = Saa'SU - ns^Ci - 

il')H{\f\)P{il). (.) stands for the ensemble average, 
and we have assumed that the frequency of the sig- 
nal is uncorrelated with respect to the angular distri- 
bution. H{\f\) and P{tl) are related to the spectrum 
f^gw(/) =dpgw(/)/Pcdln/by [p| 
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The one-sided power spectral density of the GW back- 
ground simply reads Sh{\f\) = H{\f\)/8n. P{^1) repre- 
sents the contribution to r2gw(/) from the different di- 
rections in the sky f2; in the context of our analysis, it is 
useful to expand P{Cl) in terms of the multipole moments 
Pi,n as P{n) = Y.i,rnPi^^ri^)- ^hc signal strain at 
the interferometer output is h{t) = D'^^{t) hab{x{t),t), 
where Dab = [ua{t)ub{t) — Va{t)vb{t)]/2 is the detector 
response tensor, and Ua and Va are the unit vectors point- 
ing, at each time, in the direction of the instrument arms; 
x{t) describes the trajectory of the interferometer cen- 
tre of mass, and for LISA the explicit time-dependent 
expressions can be found in The detector output 
o{t) — h{t) + n{t) is affected by the instrumental noise 
n(t), that we assume to be Gaussian and stationary, with 
one-sided noise spectral density Sn{f). 

The key idea to detect (and study the properties of) 
an anisotropic GW stochastic background is to break up 
the data set(s) of length Tobs into several data chunks of 
much smaller size r ^ T < Tobs, over which the detec- 
tor location and orientation is essentially constant, and 
construct the new signal 



S{t) = 



-\-oc 



df5,{f,t)Q{f)5l{f,t). 



(3) 



In Eq. ^ the indexes j and k label the instruments, and 
o{f;t) = J^^lj'^ dt' e~'^^^f^' o{t') is the Fourier transform 
of the data chunk of length r, centered around t. S{t) is 
the correlation function suitably weighted by the optimal 
(real) filter function Q{f) (to be determined in order to 
maximize the signal-to- noise ratio) One then looks 
for peaks in S{t) at multiples of the LISA rotation fre- 
quency l/T: because t <^ T the correlation will vary 
as LISA changes orientation. To fix ideas, we can set 
r ~ 10^ sec and Tobs ~ 10* sec. If one considers two sep- 
arated instruments, an additional time scale appears: the 



light-travel-time between the instrument centre-of-mass 
For two detectors in Heliocentric orbit at distance 
D, f ~ 500 (T'/AU) sec. As most of the correlation be- 
tween two detectors builds up for frequencies / ^ 1/^, 
there is significant signal for time-scales shorter than r. 
Since the LISA motion is periodic, one can decompose 
the mean value of S{t) as a Fourier series: 
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where each harmonics Sm reads 
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The amplitude of Sm represents the observable whose 
signal-to-noise ratio is (S/N)„j — pLm/'^m, where /im = 
\{S„i)\, and cr^ = {SmS-m) - \{Sr,i)?- It is clear that 
the noise contribution to the ensemble average (), and, 
as a consequence, the SNR of each harmonics, is different 
whether Oj and Ofc share correlated noise or not. We will 
therefore analyze the two cases in turn. 

We start by considering the LISA mission. In this 
case Oj = Ok = o. Evaluating the expectation values (.), 
and defining a suitable positive semi-definite scalar prod- 
uct [ p^ , which allows us to derive the explicit expression 
of Q{f), the SNR of the m-th harmonics becomes 
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where 
and 
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3/4 is the overlap reduction function [^,0 for 
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two co-located and co-aligned detectors, and 
5/(16^2) J^^daJdnYim{n)j:^ [Dab{a)ef{n) 
the overlap reduction functions corresponding to the Z-th 
multipole moment and the m-th harmonic. Notice that 
for co-located and co-aligned instruments 70 and 7^^^^ do 
not depend on the frequency. It is straightforward to 
show that 7;^^^ = for Z odd and I > 4. LISA is therefore 
sensitive only to quadrupole (1 = 2) and octupole (1 = 4:) 
anisotropy, and is completely "blind" to the whole set of 
odd multipole moments. This property affects both the 
signal detection and parameter estimation: LISA is not 
sensitive to the (usually strongest) dipole anisotropy, and 



2 



only two multipole moments can be exploited in order to 
extract information. 

The function j''^\f) in @ takes a particularly sim- 
ple form in the two physically relevant cases: (i) If 
h{t) ^ n{t) {e.g. for backgrounds generated by galactic 
binaries for / < 3 mHz |@) then j'^°\f) l/ll; (ii) 
For h{t) <^ n{t) fthe likely situation for primordial back- 
grounds 1), ~ [25„(/)/55„(/)]2. Using the two 

approximate expressions for j'^^\f), and assuming that 
the main contribution to the SNR comes from a frequency 
band A/ ~ /, over which S'„(/) and Sh{f) are roughly 
constant, Eq. can be cast, respectively, in the form 




where rj = Tobs A/; as reference values we consider A/ = 
10"'^ Hz and Tobs = 10* sec. he and /iims are the char- 
acteristic amplitude of the signal and the noise, respec- 
tively, per logarithmic frequency interval (see e.g. 
Notice that when he <^ /irms, the SNR is reduced by 
the factor [/ic(/)/^rms(/)]^ ^ 1 with respect to the case 
where the (isotropic component of the) signal dominates 
the noise. Eq. (^ represents the sensitivity limit of LISA 
to anisotropic backgrounds: when h^f) ^ h^ms{f)i the 
reduction of the instrumental noise does not improve the 
SNR. 

In the LISA frequency band we expect a strong 
stochastic background generated by the galactic popu- 
lation of close binary systems, mainly white dwarf (WD) 
binaries [^-^ . WD binary systems are mainly located in 
the galactic disk, and the background that they generate 
is anisotropic for an observer in the solar system. In or- 
der to gain insight into the LISA ability of detecting such 
signature we model the WD distribution as being uni- 
form within a cylinder of radius vq and height 2zo (with 
ro > zq); typical values are vq « 5.5 kpc and zq w 0.3 
kpc [||. By making this approximation the computation 

of Pim and Qm'^ simplifies considerably. For Sh{f) and 
Snif) we use the estimate provided in and ||l|, re- 
spectively, whose analytical approximations are given by 
Eq. (3.1) and (2.6) of §. Eq. (|) yields: 

This result clearly shows that LISA can unambigu- 
ously identify the anisotropic component of the stochas- 
tic background generated by galactic WD binaries. By 
exploiting this signature one could also disentangle the 
galactic contribution from the extra-galactic one, which 
is isotropic to a large extent, being dominated by sources 
at 2; 1 HJ. The large SNR at which detection can be 



made will enable us to extract valuable information about 
the distribution of binary systems in the galaxy. How- 
ever only the harmonics with to = 1 and 2 will clearly 
stand above the noise; for m > 2 the SNR is negligible. 
We have checked that the former results are weakly sensi- 
tive on the approximations that we have made: Eq. ( |ll| ) 
agrees within w 15% with the value obtained by the nu- 
merical integration of Eq. (^) using Eq. (^ over the entire 
sensitivity window and considering a more sophisticated 
source distribution. The galactic background from WD- 
WD binaries satisfies the condition hc{f) ^ hrmsif) over 
the sensitivity window where most of the SNR is accu- 
mulated. If other galactic populations in the disk gener- 
ate a stochastic background buried into the noise, LISA 
could in principle detect the signal, but at a SNR smaller 
than ( prj ) by the factor ~ [/ic(/)/^rms(/)]^, which makes 
rather unlikely to discriminate the two populations. 

LISA is very effective in detecting the anisotropy of sig- 
nals from sources located in the disk; one might wonder 
whether populations characterized by a spherical distri- 
bution (such as MACHO binary systems in the Galac- 
tic halo ||l6| ) would also produce a detectable anisotropy. 
We assume for sake of simplicity that the sources are uni- 
formly distributed within a sphere of radius re- Depend- 
ing on the value of Tc, whether it is smaller or greater than 
TGC — 8.5 kpc, the distance of the solar system from the 
Galactic Centre, the sources produce a different degree 
of anisotropy. We consider first the case of a bulge dis- 
tribution, where Tc < tqc- Assuming /ic(/) ^ /irms(/), 
and Tc — 1 kpc, Eq. (^) yields: 




We have checked that if Tc = 0.5 kpc, the SNR is es- 
sentially unaffected, and for Tc — 7 kpc the SNR be- 
comes ~ 78 and 7 for m = 1 and 2, respectively. If 
he ^ /iims the values in Eq. (^ are reduced by the fac- 
tor w 0.3 [/ic(/)/^rms(/)]^- For sources located in an ex- 
tended halo, with > tqc j the multipole moments scale 
as Pirn ^ {rGc/rcY, and Eq. (||) yields: 

(I),-°(^)"(t^)"^ <-) 

LISA can therefore provide information on a background 
generated by sources characterized by a spherical distri- 
bution. Eqs. (0), (|l|) and (|l|) show that LISA can 
indeed offer a radically new mean to explore the distri- 
bution of compact objects in the galaxy. However, the 
task of pin-pointing the exact structure of the population 
will be highly non-trivial fl^ . 

Although the former results are encouraging, they also 
suggest that it is unlikely to be able to gain significant in- 
formation on the much smaller degree of anisotropy pro- 
duced by extra-galactic generated backgrounds, for which 
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6m^ <C 1, even if he ^ /irms, see Eq. (p|). A fortiori the the intrinsic noise and the astrophysical background (i.e. 



study of the anisotropy of the primordial GW background 
is essentially ruled out. In principle one could attempt 
to detect the dipolc anisotropy induced by the motion 
of our local frame with v/c ^ 10^^ with respect to the 
Hubble flow. However LISA is not sensitive to Z = 1, and 
one would have to rely on the next "visible" multipole 
moment (1 = 2). For / = 2 we have ~ 10~^, and it 
would require a ridiculous time of integration to achieve 
SNR > 1. 

The tremendous sensitivity ensured by the space tech- 
nology, and the clear advantage of flying two indepen- 
dent LISA instruments ^ - possibly with the option 
of achieving optimal sensitivity in the frequency band 
0.1 Hz ^ f ^ 1 Hz, which is likely free from astro- 
physically generated backgrounds ~ is stimulating further 
work regarding the new range of opportunities provided 
by cross-correlation experiments [|l^,|l^. The main ad- 
vantage of this configuration for the study of the back- 
ground anisotropy is the full sensitivity to the entire set 
of multipole moments /: in general jimif) ^ Wl [ p^ . 
The value and frequency dependence of jimif) depend on 
the actual separation and relative orientation of the two 
instruments, and one could envisage to "tune" the con- 
figuration to some particular anisotropy structure. For 
two identical interferometers located at a distance D, the 
SNR reads (cfr. also |l|) 



vTobs 



dfjimiif) 



1/2 



(14) 



where J{f) and Qm{f) are the analogous of j'''^\f) and 
obtained by replacing 70 and 7^^^ with 7(7) and 
limif), respectively. 

We have estimated the SNR that one could achieve 
by cross-correlating the data sets of two identical LISAs, 
placed at a distance D = 1 AU, and for a few "random" 
orientations of the instruments. The key difference with 
respect to the LISA mission is that the odd multipole 
moments, in particular I = 1 and I = 3, become visible; 
as a consequence the harmonics with to = 1 will become 
stronger, and also the harmonic with to = 3 could be 
observable. As an example, for a disk distribution, one 
would achieve (S/N)3 ^3; for bulge and extended halo 
distributions (S/N)i would be a factor w 2 and 10 higher 
than (^2|) and (|l3|) , respectively. This will impact on the 
ability of the instrument of gaining insight on the source 
distributions, and disentangling the contributions from 
the disk, the bulge and the halo. However this particular 
aspects deserve further investigation. 

Cross-correlations would also make marginally possi- 
ble to set limits on the anisotropy of the GW primor- 
dial background, thanks to the sensitivity to the dipole 
moment. Let us consider a primordial background with 
P(f2) = 1 -f -^0 + 0(/3|j), where the parameter Pd 
measures the degree of anisotropy and the unit vector 
rio defines the direction of maximum anisotropy. As- 
suming that the primordial spectrum overwhelms both 



iJqq rjgw > 10 ") the parameter must satisfy Pd > 
8 X 10-3 (robs/10^sec)-i/2 (A//10-3Hz)-i/2 ^^^^^j. 
to detect the dipole component. The dipole anisotropy 
induced by the motion of the local reference frame with 
respect to the Hubble flow corresponds to Pd ^ 1Q~^; 
therefore only an instrument with optimal sensitivity 
in the band / ^ 0.1 Hz would be able to pick it up. 
This band presents also the considerable advantage of 
being completely transparent to primordial GW back- 
grounds 1^. Assuming that one can perfectly remove 
the induced dipole anisotropy from the data using inde- 
pendent information provided by CMB experiments, one 
could then be able to test an intrinsic statistical depar- 
ture from isotropy compatible with the measurement of 
the quadrupole anisotropy of the CMB temperature Ml . 
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